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Abstract 



> 

in 

f~^ . Invariants at arbitrary and fixed energy (strongly and weakly conserved quantities) for 2-dimensional 

CNJ ' Hamiltonian systems are treated in a unified way. This is achieved by utilizing the Jacobi metric ge- 

_|I , ometrization of the dynamics. Using Killing tensors we obtain an integrability condition for quadratic 

("■^ ' invariants which involves an arbitrary analytic function S{z). For invariants at arbitrary energy the 

QQ ' function S{z) is a second degree polynomial with real second derivative. The integrability condition then 

f^ , reduces to Darboux's condition for quadratic invariants at arbitrary energy. The four types of classical 

quadratic invariants for positive definite 2-dimensional Hamiltonians are shown to correspond to certain 

conformal transformations. We derive the explicit relation between invariants in the physical and Jacobi 

^^ ■ time gauges. In this way knowledge about the invariant in the physical time gauge enables one to directly 

$_( ' write down the components of the corresponding Killing tensor for the Jacobi metric. We also discuss 

5^ , the possibility of searching for linear and quadratic invariants at fixed energy and its connection to the 

problem of the third integral in galactic dynamics. In our approach linear and quadratic invariants at 

fixed energy can be found by solving a linear ordinary differential equation of the first or second degree 

respectively. 
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1 Introduction 

The quest for integrable systems is still one of the main areas of interest in classical dynamics both per 
se and for the applications in the related fields of celestial mechanics, accelerator physics and so on. In 
galactic dynamics a long-standing and to this day unresolved issue concerns the celebrated third integral. 
Briefly, the observed motion of stars in the galaxy has a certain regularity which indicates the existence of 
a third invariant in addition to the energy and angular momentum. This picture is confirmed by numerical 
simulations. In modern terminology, the motion is non-chaotic contrary to what one would expect for a 
generic potential. A particularly striking numerical experiment was made by Henon and Heiles [1] in 1964 
using a model potential now known as the Henon-Heiles potential. They showed that the motion in that 
potential is regular with invariant tori up to a certain energy value Eq. At this energy the invariant tori 
abrubtly begin to break up and the motion becomes chaotic in increasingly large parts of the phase space. 
One conclusion is that there can be no invariant which commutes with the Hamiltonian. On the other hand 
the very high degree of regularity at low energies, E < Eq, seems to indicate the existence of an additional 
invariant at those energies. Such an invariant, if it exists, must necessarily depend on the energy with some 
kind of singular behaviour at E = Eq. However, no such invariant is known, for any value of the energy. 
One of the purposes of this paper is to discuss the possibility to find invariants at fixed energy values. 

An essential aspect of integrability is therefore to check for the existence of integrals of the motion in 
addition to the energy, namely phase-space functions invariant with respect to the phase flow. The search 
for additional invariants has therefore been actively pursued for more than a century (Bertrand, 1852 |2]i 
Darboux, 1901 [3]; for a complete review of the matter see Hietarinta [4]). A less known aspect of the question 
is the relation between the invariants that most frequently appear in the applications, namely polynomials in 
the momenta, and Killing tensors of the configurational manifold of the dynamical system. This relation is 
in the form of a straightforward correspondence if the Jacobi approach to the geometrization of the dynamics 
of a conservative system is adopted. 

The framework of the Jacobi geometrical formulation of the dynamics offers several conceptual and 
technical aids that give the possibility of shedding new light on a number of aspects of integrability. From the 
technical point of view, it allows the use of the powerful tools of Riemannian geometry on the configurational 
manifold of the corresponding natural Jacobi formulation of the dynamical system. From the conceptual 
point of view, it treats invariants at fixed energy (the so called configurational invariants of Hall [5]) on the 
same footing as the invariants at arbitrary energy, thus realizing a unified geometric approach. 

One aim of the present paper is to demonstrate the fruitfulness of the geometric approach by providing 
the generalization of Darboux's conditions to include the configurational invariants. In the particular but 
fundamental case of two degrees of freedom, the Killing equations for second rank Killing tensors can be solved 
in full generality, resulting in all the known cases of integrability at arbitrary energy, both the classical ones 
already obtained by Darboux [3] and those recently found to exhaust all the possibilities in two dimensions 
(Dorizzi et al. [6]). 

The substantially larger family of solutions of the Killing equation in the fixed energy case allows one 
to find large new classes of constrained integrable systems, with a backwards procedure going from the 
knowledge of the form of the invariants to the assessment of the structure of the potentials admitting 
constrained integrability. Applying a new prescription for transforming Hamiltonian symmetries between 
different time gauges, the relation between the invariants in the physical and Jacobi time gauges is explicitly 
determined. The relation between the invariants also provides a shortcut to calculate the Killing tensor from 
the knowledge of the invariant in the physical time gauge. 

2 Geometric representation of the dynamics 

The geometrization of the dynamics in terms of the Jacobi geometry has been known for about a century 
and can be found in some text books [3 [HI [9] . In spite of this it is not widely known or used and therefore 
we outline the main ideas in this section to make the paper self-contained. Most Hamiltonian systems 
of physical interest have a geometric kinetic energy part, that is to say that the kinetic energy is a non- 



degenerate quadratic form in the momenta p^ 

T = ^/i"Vp/3 , (1) 

where /ia/3 is a function of the configuration variables g" . The Hamiltonian itself then has the form 

H = T + V{q)=E . (2) 

The independent variable t is often but not always the time. For simplicity we shall refer to the independent 
variable as the time in this paper. For any given energy E of the system we can use the Hamiltonian 

H = H -E , (3) 

to represent the dynamics provided that we impose the constraint 

H = . (4) 

For any such zero energy Hamiltonian we can reparametrize the system by introducing a new time variable 
tj\f defined by the relation 

dt = JV{p, q)dt^r , (5) 

together with a redefined Hamiltonian 

n^f^^f{p,q)H^^fT+^f{v-E)^o. (6) 

The new Hamiltonian will then give the same equations of motion on the constraint surface Ti^v — (see 
e.g. [in])- We shall use the term lapse function for Af{p,q) which defines the independent variable gauge. 
This usage is borrowed from applications in general relativity where the lapse gives the rate of physical time 
change relative to coordinate time (see e.g. jllj ). The lapse function can be taken as any non-zero function 
on the phase space. 

The dynamics of the system can be represented in a purely geometric formulation by exploiting the 
reparametrization freedom. The passage to the geometric representation is accomplished by defining a new 
time variable tj [Jacohi time) by the lapse choice 

M^Mj^[2{E~V)]-' . (7) 

Note that for a positive definite kinetic energy E — V is always nonnegative in the physically allowed 
region. The corresponding Hamiltonian has the form Tij^f = [2(i? — V)\^^T — i. Thus TLj^ has a constant 
potential energy which can be subtracted without affecting the equations of motion. This leads to the Jacohi 
Hamiltonian^ Hj, defined by 

Hj = \Njh^"p^P[3 = \r^Pc.Pp , (8) 

where J"'^ is the Jacohi metric of the system. This defines a geometry which contains all information about 
the dynamics. The dynamical metric is given in covariant form by 

J^p = 2{E - V)h^0 . (9) 

The Jacobi metric is therefore conformally related to the original kinetic metric hnp. It follows that the 
dynamics of the system ^ is equivalent to geodesic motion in the Jacobi geometrjo 

ds^ ^ Japdq"dq>^ . (10) 



^ It may happen that there are other geometries which can also be used to represent the dynamics. In that case the dynamical 
geometry defined according to the prescription given here is not unique. For example an inequivalent dynamical geometry may 
sometimes be obtained by performing a suitable canonical transformation in some time gauge (see e.g. I12p . 



Note that the Jacobi geometry depends on the energy parameter E. This means that in general the Jacobi 
geometries corresponding to different energy surfaces are inequivalent. The complete geometric representa- 
tion of the dynamics of the system is therefore given by the geodesies of a l-parameter family of geometries. 
Because of the form of the conformal factor in the Jacobi geometry ^ the geometric representation is 
only valid locally in configuration space at points where V y^ E. In general the equation V = E defines 
a non-empty energy surface in configuration space for a given value of E. Exceptions occur if e.g. T^ > 
throughout the configuration space and E <0. In this paper we are only concernced with the local existence 
of invariants. Therefore the failure of the Jacobi geometry to represent the global dynamics does not affect 
our analysis. 

3 Positive definite two-dimensional Hamiltonian systems 

To study 2-dimcnsional systems it is very helpful to use variables which are null (lightlike) with respect to the 
dynamical metric. Such variables are naturally adapted to the action of the conformal group which plays an 
essential role for 2-dimcnsional systems. The indefinite (Lorentzian) signature case was discussed by Rosquist 
and Uggla [13j . The null variables for that case are real and the conformal group can be parametrized by 
two arbitrary real functions of one variable. Our approach is to use the method of [13 to treat the case of 
a positive definite dynamical metric. 

We consider a general 2-dimensional positive definite dynamical metric written in the manifestly confor- 
mally flat form 

ds^ = 2G{x,y)idx^+dy^) . (11) 

We seek a condition on G which guarantees the existence of a second rank Killing tensor. To that end we 
follow [13] as closely as possible and introduce null variables which in the positive definite case are complex 

z ~ X ~ iy . 

Here and in the following a bar is used to denote complex conjugation. The metric then becomes 

ds^ = 2G{z,z)dzdz . (13) 

Note that although our notation is similar to that in pO?] the variables z and z are complex in the positive 
definite case as opposed to the Lorentzian case where they are real. 

When doing calculations it is convenient to employ either an orthonormal frame ui^ using hatted upper 
case Latin indices (/, J, K, . . . — 1,2) in terms of which the metric is written as 

ds^ = {Jf + {Jf , (14) 

or a complex null frame fl^ using upper case Latin indices {I, J,K, . . . — 1,2) with the metric written as 

ds^ = 2f7irj2 . (15) 

The frame components are given by 

J = {2Gy^^dx , n^ = G'/^dz , 

. (16) 

u^^{2Gy/^dy, n'^G^'^dz, 

The two frames are related by i7^ — e^ fto^ or lo^ — ej^il^ where e/^ is the transposed matrix inverse of e^ j. 
The transformation matrix is given by 



(e^-) = f ^ , (e/') = f f ] . (17) 
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The simplest symmetry of the geodesic equations is known as a Killing vector field (see e.g. [II]). A 
second rank Killing tensor is a symmetric tensor Kij — K(^jj^ satisfying the equation 

KiIJ;K) = . (18) 

The existence of such a tensor is equivalent to the existence of a polynomial invariant of the second degree 
K^'^PjPj for the geodesic equations. Killing tensors and second degree invariants are the natural general- 
izations of Killing vectors K' and the corresponding first degree invariants K^pj. Any Killing vector gives 
rise to a Killing tensor K/jKjy Such a Killing tensor is said to be reducible. An important property of 
the second rank Killing tensor equations is that they can be decomposed in conformal (traceless) and trace 
parts according to 

P{IJ;K) - 2^(IjP K);L = , 

K.J = ~P-\,j , 
where the Killing tensor itself is decomposed in a conformal part P/j and the trace K = K^ i according to 

Kij - Pij + \Khij . (20) 

Referring to the vector P/ = P^ i-^j as the conformal current it follows from (|19p that for any given conformal 
Killing tensor the equation for the trace can be solved if the integrability condition 

P[i,J] = (21) 

is satisfied. The procedure to solve the Killing tensor equations is therefore to first solve the conformal 
Killing tensor equations and then check if the integrability condition (PT|) for the trace can be satisfied. In 
fact it turns out that the conformal equations can easily be solved leaving the integrability condition as the 
remaining equation to study. 

Let Kjj be the components of a Killing tensor with respect to the orthonormal frame w^. Then the null 
frame components are given by Kij — ej-'^'^ ej^ K j^j^j which gives 

i^n-^(i^ii-i^22)-*^i2: 

^22-^(i^ii-i^22)+»^i2: (22) 

where K — K^ j is the trace. Since the orthonormal components Kjj are by definition real it follows that 
K22 — Kii. Therefore we can represent the null frame components of the conformal part of the Killing 
tensor by a single complex function S'(z, z) according to (cf. [13j ) 

the trace K being a real function. The conformal Killing tensor equations are formally identical to those 
obtained in the Lorentzian case |13| 

Cm = G^'^S^, = , C222 = Gi/25 , = . (24) 

It follows that Pa/at is a conformal Killing tensor precisely if 5 is a function of z only. One way of stating this 
result is that the equations ([M)) for the conformal Killing tensor coincide with Cauchy-Riemann's equations 
for certain rescaled linear combinations of the tracefree parts of the Killing tensor. This result was implicit 
in [TS]. Thus the conformal Killing tensor has the form 



The trace equations can then be written in the form 

if,, = -2S{z)G^, ~ GS'iz) , 

K^g = -2S{z)G^, - GS'iz) , ^ ' 

the derivative K^g being given by taking the complex conjugate of the above equation. 

As in the Lorentzian case there is an infinite-dimensional family of conformal Killing tensors in this case 
parametrized by the single complex analytic function S{z). Following the procedure in [13j we now write 
down the integrability condition for the conformal current. This is exactly what we need to guarantee the 
existence of a Killing tensor. In the positive definite case the integrability condition becomes 



PlU2] = - G-'G,,,Siz) + G-'G,ggS{z) - ^G-^G,.S'{z) + '^G-'G,gS'{z) 
-ls"iz) + ^S"iz)^0, 



(27) 



or in a more compact form 

P[i;2] = -i 3{2G"'G,..5(z) + 3G-^G,,S'iz) + 5"(z)} = , (28) 

where 9 denotes the imaginary part. As in J13j we use a conformal transformation to standardize the 
frame and coordinate representation of the conformal Killing tensor. To that end we introduce the new 
null frame il^ — Bfl^, (l'^ = B~^^^ and a new complex coordinate by the transformation w = H{z) with 
inverse z = F{'w). By choosing B = {F' {w) / F' {w))^/"^ we ensure that the new frame ^^ has the coordinate 
representation (l^ = G^^^dw, Cl^ = G^/^dw where G = \F'{'w)\'^G is the new metric conformal factor, 
ds^ — 2Gdwdw. It follows that the conformal Killing tensor components in the new frame are given by 

P„ = B-'Pn = [H'{z)]^Si-z)G , 
P22^B^P22^[H'iz)]^Siz)G. 

The coordinates are standardized by choosing a conformal transformation function H{z) which satisfies 
[H'{z)]~^ = S{z) implying that Pn = P22 = G- Note that det(P/j) is always positive. Therefore there is 
only one type of Killing tensors for positive definite Hamiltonians unlike the situation in the Lorentzian case 
(see [13] )■ For a standardized conformal Killing tensor the integrability ^7\ condition simplifies to 

Writing w = X + iY the solution of this equation is 

G = QiiX) + Q2iY) . (31) 

This is the usual form of the potential in separable coordinates. 
The trace equation [see (fT9)) ] can then be written 

""'" = "'^•" ' (32) 

K^Y = 2G,y , 

with the solution 

K^-2Qi{X) + 2Q2{Y). (33) 

The integrability condition (j27p can be interpreted as follows. Given any analytic function S{z) there is a 
family of potentials (the solutions of ([27| ) which is integrable at zero energy. For future reference we note 
that the original potential is given in terms of the standardized potential by the relation 

G^\S{z)\~^G. (34) 



We wish to express the invariant Ij = K^'^^pmPn in terms of the coordinate momentum components. 
The subscript J is used to distinguish the invariant in the Jacobi time gauge from the invariant / in the 
physical time gauge (see section [B]). By (|16l) the nuh frame components are related to the complex coordinate 
components by (pi,p2) = G^^''^{pz,Pz)- The invariant can then be written in complex coordinates as 

/,/ = 5(z)p,2 + siz)p,^ + G-^KpzPz . (35) 

Using the relation 

Pz = -^[Px +iPy) , 

implied by ([T2|) we obtain the invariant in real coordinates as 



(36) 



I J - I'iiiSKPx^ - Py^) + ^iS)p.Py + ^G-'KipJ+p,/) . (37) 

Note that the last term is equal to KHj. 

4 Arbitrary energy invariants 

A common situation is that one is interested in invariants which are valid for arbitrary values of the energy. 
Invariants of that type arise if the integrability condition itself is independent of the energy. This happens 
precisely if 3{S"'(0)} = 0. Since S"{z) is an analytic function it follows that S"{z) is a real constant and 
therefore S{z) must be a second degree polynomial 

S'(z) = M/^(z) := az^ + /3z + 7 , (38) 

where a is a real constant and j3 and 7 are complex constants. The integrability condition (|27p then simplifies 
to 

'^{2G^zzS{z) + iG^zS'{z)}^Q . (39) 

This is nothing but Darboux's condition for quadratic constants of the motion [3l [4j written in complex 
variables. The general integrability condition (|27p therefore generalizes Darboux's condition to include also 
quadratic invariants at fixed energy. 

For a given analytic function S{z) the potential may be written as G{X, Y) = \F'{'w)\^'^{Qi{X) + Q2{Y)). 
In the arbitrary energy case we have S{z) = W{z) = [H'{z)\~'^. The equation H'{z) = [iy(z)]~^/'^ can then 
be integrated resulting in (without loss of generality we may put a = 1 if a 7^ since p9|) is invariant under 
a real scaling of S) 



log(VT4^ + Z + (5/2) + const , (a = 1) , 
H{z) = { 2/3-V/':^2 + 7 + const , (a == 0, /3 7^ 0) , (40) 

7^1/^2 + consi , (a = 0,/3 = 0) . 



Inverting these relations we must distinguish between the four cases 

(i) a = 1,A:= v//3V4-a7 7^0 

(ii) a = l,A = 0, 

(iii) a = 0, /? 7^ , 

(iv) a = 0, /? = , 

giving rise to the following explicit formulas for F{w) 



(41) 



F{w) 



' Acosh(w--u;o)-/3/2 , (i) 

\l3{w ~ wjf - 7//3 , (iii) 
^ 7"'"/^w — Wo . (iv) 



(42) 



where wq is an arbitrary constant. We are primarily interested in proper conformal transformation, i.e. 
S{z) y^ 1. Therefore we are free to perform translations and rotations in z and w to simplify formulas. In 
particular we may put wq equal to zero by a translation of the origin in the w-plane. We can also rotate the 
2;-plane to transform A to a positive real number in case (i). Likewise (3 may be assumed real in case (iii) 
and 7 can be taken as real and positive in case (iv). The argument we used to set a = 1 can then be used 
to set /3 = 4 and j = 1 for convenience in cases (iii) and (iv) respectively. We also translate the z-origin 
obtaining finally 

{A cosh w , (i) 
^^2 ' J--\ (43) 

w , (ni) ^ ^ 

w , (iv) 

where A is now to be understood as a positive real number. Differentiation of the expressions (|43p yields 




^'H= L' n (44) 



implying that the conformal transformation factors take the forms 

A2 (sinh^ X + sin^ Y) = ^{r^ + A^y - 4A^x^ , (i) 

iswiHF'wPH:;^r;;^,,,,, <;;;, (45) 

1 , (iv) 



where r := a/o^+j/^. 

In case (i) the conformal coordinate transformation is given by 

a; = A cosh X cos Y , 
y — A sinh X sin Y . 



(46) 



This is the transformation used to define the Stackel potential [16] which was also discussed by Darboux [3]. 
Case (i) therefore gives Stackel's classical integrable potential. From ((45|) and (p4|) we find that the potential 
is given by the well-known formula [4] 

^r , Qi{X{x,y)) + Q2{Y{x,y)) _ A{^{x,y)) + B{r^{x,y)) 

^""'^^ A^[smh^X{x,y) + sm^ Y{x,y)] Cix,y) + Tjix,y) 

where 

^{x, y) := 2A2 sinh^ X = r^ ~ A^ + ^{r"^ + A^^f - AA^^x"^ , 

r/(a;, y) := 2A2 sin^ Y = -r'^ + A^ + ^(r^ + A2)2 - 4A2a;2 . 

and the functions A and B are arbitrary functions of their arguments. The case (ii) conformal transformation 

is given by 

X = e-^ cos Y , 

X ■ ^ (49) 

y = e sm Y . 

In particular it follows that Y — a.i:cta.n{y / x) ~: is a polar angle. Using ([^5]) it then follows that the 
potential can be written as 

Gix,y)=A{r)+r-^B{(t)) , (50) 

where A and B are arbitrary functions. This is again a classical integrable case [J known as the Eddington 
potential. Case (iii) is characterized by the conformal transformation 

x^X^ -Y^ , 

51 
y = 2XY . ^ ' 

8 



conformal transformation conformal Killing tensor component „ . ,. 

case ,,/ ^ „, s rrr// M-9 Separating coordmatcs 

w~H[z) b(z) — [H(z)\ 

~(i) lnz±Vz^-A^ z2- A2 Elliptical 

(ii) Inz z^ Spherical 

(iii) ^/z Az Parabolic 

Table 1: Conformal transformation functions giving rise to systems with a second invariant 
which is linear or quadratic in the momenta. 



Using (|45p this gives immediately the likewise well-known classical integrable potential [3] 

G{x,y) = [A{r + x)^B{r~x)]/r , (52) 

where again A and B are arbitrary functions. Finally in case (iv) the conformal transformation factor is 
unity and so the potential is simply given by the explicitly separated form 

G(x,v) = A{x)+B{y) , (53) 

in terms of the arbitrary functions A and B. This completes the list of the four classical cases. The three 
nontrivial conformal transformations giving rise to systems with a second linear or quadratic invariant are 
tabulated in table [1] 

5 The Killing vector subcase 

Although Killing vectors correspond to reducible second rank Killing tensors it is nevertheless worthwhile 
to give a separate treatment of that subcase. As will be shown below the integrability condition for Killing 
vectors is a first order differential equation. In some applications it can advantageous to investigate this 
simpler case before tackling the second rank Killing tensors. A Killing vector can be described by a function 
Z{z,z) according to (cf. [13]) 

The Killing vector equations then become 



Kil.l) = i^(2;2) - Z^, = , (55) 

i^(i;2) = \ {Z,. + Z,-z + G-^G,,Z + G-^GrzZ) = . (56) 

This shows that Z depends only on z and is therefore a complex analytic function. It follows that the 
remaining Killing vector equation (j56p reduces to a form analogous to the integrability condition for the 
second rank Killing tensors ((27)l . It can also be written in the form 

K^^.2) = G-^3i{GZ'{z) + G^,Z{z)} = G-^3l{{GZ)^,} = . (57) 

Any Killing vector gives rise to a second rank Killing tensor given by 

Kij = KiKj + chij , (58) 

where c is an arbitrary (real) constant. The conformal Killing tensor components then become 

Pn = i^ii = {K,f = GZ^ , 
P22 = i^22 = {K2f = GZ\ 



Referring to ([^5]) it follows that the conformal Killing tensor is determined by the analytic function 

S^Z^ . (60) 

As in the second rank case, invariance at arbitrary energy involves a further restriction coming from the 
requirement that (|56p should be invariant with respect to energy redefinitions. This means that we must 
impose the condition 

3fi[Z'(z)] = , (61) 

leading to the simplified Killing vector equation 

G^^Z + G.jZ^O . (62) 

where 

Z{z) = U{z):=ihz + 5 , (63) 

and h and 5 are real and complex constants respectively. It follows that the corresponding second rank 
conformal Killing tensor is determined by 

S{z) = [U{z)f = ~}?z^ + 2ib5z + 5^ . (64) 

Thus S{z) is an even square for a reducible Killing tensor so the condition A^ = /3^/4 — 07 = is always 
satisfied. 

6 Transforming Hamiltonian symmetries between different time 
gauges 

In this section we address the problem of how an invariant is affected when we transform from one time 
gauge to another. Suppose we are given a Hamiltonian H with a second invariant / so that {/, H} = 0. In 
another time gauge the Hamiltonian is given by TLj\/ = MH. = N(H — E) where M = dt/dt. Of particular 
interest for the purposes of this paper is the transformation between the physical and Jacobi time gauges. 
In that case, going from the physical time gauge to the Jacobi time gauge, dt ~ Mdtj, we have Af = \2V\^^. 
Since the Poisson bracket 

{I,H^} = H{I,Af} , (65) 

is in general non- vanishing off the zero energy surface Ti. = 0, the original invariant does not have a vanishing 
Poisson bracket with the Hamiltonian in the new time gauge. Borrowing usage from Dirac's theory of 
constrained Hamiltonians we may say that the invariant is only weakly conserved in the new time gauge. 
For a linear invariant, e.g. the Killing vector case, one can always choose variables such that I — px and x 
is a cyclic variable in the Hamiltonian. Then for the Jacobi time gauge, Af depends only on the remaining 
variable and hence {/, Hj^f} = implying that / is strongly conserved in the new time gauge. Suppose now 
that we have a 2-dimensional system with a second invariant given in the physical time gauge. Then, the 
system is integrable and we can express the lapse explicitly as a function Af — f{t) of the physical time 
t. Integrating the relation dt = dt/ f{t) then gives the new time as an explicit function of the old time at 
least up to a quadrature. In principle, the solutions can consequently always be transformed to another 
time gauge. One therefore expects that the integrability properties of a dynamical system are independent 
of the time gauge. In order to exploit fully the geometrical formulation of mechanics it is desirable to find a 
corresponding invariant which is strongly conserved in the Jacobi time gauge. As we shall see this is actually 
possible at least in the cases considered in this paper. 

To find the invariant /v in the new time gauge it is convenient to use an ansatz of the form 

Im = I + RHm . (66) 
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We look for a condition on R which guarantees that /v is a constant of the motion for Hj^. To that end we 
form the Poisson bracket 

{i^,n^} ^ Mn^[{R,n} - {i,M-^} -n^lR,^'^}]. (67) 

Requiring this expression to vanish off the zero energy surface yields the condition 

{R,n} = {i,M-^} + nM{R,M'^}- (68) 

This is a necessary and sufficient condition for /a/ to be a constant of the motion for Tij^f. In general a 
solution to this partial differential equation would be difficult to find. For our purposes, however, it turns 
out that we can actually find a solution by a simple procedure which we now outline. Suppose now that 
J is a quadratic invariant. Then if A/" is a function on the configuration space, the bracket {I,Af~^} is a 
linear function of the momenta. Suppose we find a function R on the configuration space which satisfies 
{R,Ti} = {/,A/'^^}. Then since Hm{R,-^^^} = we have a solution of ([55]) . In this way the procedure 
to find the invariant in the new time gauge is reduced to calculating the Poisson bracket {I,Af~^} and 
finding a function whose time derivative coincides with that bracket. Although we do not know under which 
conditions this procedure works it does work in the cases considered in this paper. To summarize we first 
calculate the function {I,Af~^} and check whether it can be expressed as a total time derivative of some 
function R. If in addition {R,J\f^^} — then R is the required function which satisfies 



7 The relation between the quadratic invariants in the physical 
and Jacobi time gauges 

We wish to see how the invariant p7p appears in the physical time gauge in terms of the conformal function 
S{z) and the trace K. To that end we use the procedure outlined in section [S] going "backwards" from 
the Jacobi time gauge to the physical time gauge. Referring to section [2] the starting point is now the 
Jacobi Hamiltonian H = Hj = jG~^{px^ ~^ Py^) considered at the energy value 1/2 so H. — Hj — 1/2. 
The transformation from the Jacobi time to the physical time is then given by dtj — Mdt where N = 2G. 
According to (pS)) we write the physical time invariant as / = /j + RAfTC. Using the prescription in section 
ini R should satisfy the equation {Ij,Af~^} = {R,Hj}. Calculating the left hand side {Ij,N'~^} with Ij 
given by the expression (1571) gives 

{Ij,M-'} = G-^SGG,, + ^KGrz)Pz + G-\SGG,, + \kG,,)p, . (69) 

Assuming that i? is a function on the configuration space the right hand side is given by 

{R, Hj} = G-HR,-zPz + R,zPz) ■ (70) 

Comparing equations ([69]) and ([70]) then leads to 

R ^ = SG-^G 2 + -KG-^G , , 
i?,2 = SG-^G^, + -KG-^G^, . 

At this point it is convenient to introduce a function Q — R+ {1/2)KG^^. The equations ([7T|) then reduce 
to 

Q.z = ~7:^,z , 

2 (72) 
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where we have used the trace equation (|26p . The integrabihty condition for this equation coincides with the 
condition for integrabihty at arbitrary energy, QS"{z) = 0. Using (|38|) we then find that the solution of ((72)) 
is given by 

Q = -azz- -/3z- -l3z + Qo , (73) 

where Qo is a (real) integration constant. 

Collecting our results we find that the expression for the invariant in the physical time gauge is 

I = Iconf + \q{p.^ + Py^) + ^-GQ , (74) 

where 

Iconf = P'''''PMPN = ^3fi(5)(p,2 _ py-2) + ^{S)p.,py , (75) 

is the conformal part of the invariant. It follows that this part is the same in both time gauges. It is also 
seen that Q can be identified with the trace of the physical invariant with respect to the physical metric Sab 
where a and b are coordinate indices taking the values x and y. The relation between the invariants provides 
a shortcut to calculate the Killing tensor from knowledge of the physical invariant. We shall outline this 
procedure and then illustrate with an example. Let the physical invariant be given by an expression of the 
form 

I^Q'''{x,y)paPb + f{x,y) , (76) 

The physical Hamiltonian is given by 

H=^6-'paPb + Vix,y), (77) 

We can read off the function Q by Q — 6"'^Qab- Taking the conformal part and identifying with ([75]) we 
obtain 

S = 2(P^^ + iP'^y) , (78) 

where P"'' = Q"'' - \Q5''^ and P^/^ = -P^^. Finally we obtain the KiUing tensor trace as K ^ 2(/ + GQ) 
putting G = E -V. 

Let us now illustrate the above results by an example. We take the Kepler potential V — — /i/r where 
r = va;^ + y^ with its well-known non-trivial quadratic invariant, the Laplace-Runge-Lenz vector (see e.g. 
[TT] ) with components 

1 2 \ X 
Ll ■■= e^r = -[Xpy ~ yPxPy) , 

1 I '™' 

L2 := ey = -{ypx'^ - xpxPy) . 

/i r 

The homogeneous and inhomogeneous parts of the invariants Li and P2 are consequently given by 

and 

nab _ ( y/fJ- -x/{2^i)\ _ I , , 

^m-\-^l(2y) ) ' /(2)--yA- (81) 

It follows that the conformal Killing tensor components are given by 



^m = -^m = -^/(2m) , P,^i^ = -y/(2A.) , ^ ^ 

(82) 

(2) = -^(2) "" ^ 



Pf.^. - -PJl - 2//(2/i) , P^) = --^/(Sm) , 



while the traces are given by 

if(i) = 2(Elii)x , X(2) = %Elii)y . (83) 
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8 Applications to integrability at fixed energy 

An important lesson to be learned from the present work is that integrability at fixed energy and arbitrary 
energy (weak and strong conservation laws) are just two aspects of the same phenomenon. In particular a 
fixed energy invariant in the physical time gauge corresponds to an arbitrary energy invariant if the system 
is geometrized by going to the Jacobi time gauge. To illustrate how this works in practice we give a few 
examples in section [8 . II of fixed energy invariants beginning with the Kepler potential. Surprisingly, we find 
two apparently unknown linear invariants at zero energy for the Kepler problem. It is remarkable that it 
is still possible to discover new properties of such a simple and well-known system. This is in fact a sign 
of the power of the geometric formulation of dynamical systems. In section 18.21 we discuss how conformal 
transformations can be used to generate systems which are integrable at fixed energy. 

8.1 Degeneracy of the Laplace-Runge-Lenz vector at zero energy 

Consider now the Kepler potential, V — —fJ-/r, using polar coordinates defined hy x — r cos (p, y = r sin (p. 
In this case we have one linear invariant at arbitrary energy, the angular momentum p^. We are interested 
in finding out whether there exists another linear invariant at some fixed energy value. If this is the case 
the Jacobi metric ^ has two Killing vectors. However, a 2-dimensional space with two Killing vectors must 
necessarily also have a third Killing vector ([T51, Theorem 8.15) and the geometry is then a space of constant 
curvature. This can easily be determined by computing the scalar curvature 

(2)^ = 2G-^G,,i - 2G-^G,,G,i , (84) 

of the metric p^ and checking if it is constant. For the Kepler potential we have G = E + ^{zz)~'^/'^ and 
the scalar curvature becomes 

(2)d_ -^^ CS^;^ 

2[^i + E{zzY/^f ■ ^ ' 

This shows that the Jacobi geometry has constant curvature only if _E = and then the geometry is actually 
flat. Of course the flatness of the Jacobi geometry in this case also follows directly from the form of the metric 
since G is then a product of functions of z and z. The two extra Killing vector flelds can immediately be 
written down if we introduce Cartesian coordinates, (X, Y), for the Jacobi geometry by the transformation 
w = X + iY = \/2z leading to the manifestly Euclidean form 

ds/ = 4^(dX2 + dY'^) . (86) 

The relation to the original coordinates is given by the parabolic transformation 

x^\{X^-Y^), y^XY. (87) 

This transformation does not have a unique inverse. However, in the region X — ^{w) > we may select 
the inverse transformation to be 



X = \/r + X , Y — (sgny)-\/?' — x . 

It follows that the Kilfing vector flelds are K(i) = d/dX, Kf^2) ^ d/dY and 1^(3) = ~Yd/dX + Xd/dY = 
2d/d(p. The two extra Killing vector flelds are thus the translation symmetries iCji) and i^(2)- The corre- 
sponding invariants are 

h ■■= Px = ir + xfl'^v. + (sgnz/)(r - xfl'^Vy = (Zr)^/^ cos(0/2)p, - (2/r)i/2 sin(0/2)p^ , 
h ■■= PY = -isgny){r - xf'^Vx + (r + xfl^Vy = i'^rfl^ sin((/./2)p, + (2/r)i/2 cos(0/2)p0 . 

These two invariants are related by the quadratic formula I-^ + 1-2^ = 8/ii?j where H'^j = —{2V)^^T = ^ is 
the Jacobi Hamiltonian for the zero energy system. In terms of the physical Hamiltonian, the corresponding 
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relation is /i^ + /2^ = 4/i(l — H/V) = 4/^. In fact here we have the key to the physical interpretation of /i 
and l2- To see this let us introduce an invariant (Jjq at zero energy by the relations 

h = -2fi^/'^sin{(l3o/2) , I2 = 2^i^^^cos{(j)o/2) . (90) 



Now using ([89]) to solve for the radial momentum yields Pr — •\/2/i/r sin[((/) — (J)q)/2]. Inserting this value 
into the Hamiltonian constraint H — and solving for r gives the familiar relation 

r = EHI^ (91) 

1 + cos{(t) - 0o) 

This shows that (f>o is nothing but the angular integration constant. It follows that changing the value of /i 
(or I2) only affects the parametrization of the orbit while leaving the orbit itself invariant. From this point 
of view these invariants are gauge symmetries of the Kepler system at zero energy. 

The invariants /i and I2 are in fact closely related to the components of the Laplace-Runge-Lenz vector. 
Expressing those components as 

Bx = ecos(/>o , By = esin(/)o , (92) 

where e is the eccentricity (e = 1 at zero energy) and comparing with ()90|1 it is evident that 

L, = -(2a.)-1/i2 + 1 , L2 = -(2m)-1/i/2 . (93) 

This can also be seen directly by calculating for example /i^ from the expression given in ((89|) with result 

h^- = 2{r + x)H + 2fi{l-Li) , (94) 

where we have used ([7^ . Now solving for Li at zero energy gives again the first of the relations ([M|) . 

The relations ([M]) imply that the second rank Killing tensors corresponding to the components of the 
Laplace-Runge-Lenz vector are reducible at zero energy. The invariants in the Jacobi time gauge can be 
found from the relations ([82]), dHS]), dZH]) and ([371). This gives 

Ji = ~{2fiy^x{p^'^ - p,/) ~ fi'^yp^py , 

J2 = (2/i) y{p^ ~Py ) - /I Xp^Py , 

where Ji and J2 are the Jacobi invariants corresponding to Li and L2 respectively. From ([M| we then find 

the relations 

.h^~{2fi)-'h' + 2Hj , 

, (96) 

J2 = -(2/i)-i/i/2 . 

The reducibility of the Killing tensors KM-^ and KMJ^ corresponding to Ji and J2 is therefore expressed by 



the formulas 






MN 

(97) 



We also wish to understand the commutation relations for the fixed energy invariants in the physical 
time gauge. To facilitate the calculations the physical Hamiltonian is first expressed in terms of the Jacobi 
Hamiltonian by 

H = T + V -E={V -E){l-2Hj) , (98) 

where we have used T — 2{E — V)Hj. We can now exploit the fact that the invariant commutes with the 
Jacobi Hamiltonian to obtain 

{/, H} = {/, H} = (1 - 2Hj){I, V} = n{V - Ey'{I, V} . (99) 

This relation shows that we only need to compute the Poisson bracket with the potential. It also follows 
that the bracket {I, H} in general depends linearly on Ti,. 

Using ([M)) to calculate the brackets for the Kepler invariants we find 

{Ii,H} =r-^{r + xfl'^H = \/2r"^/^ cos(</>/2)iJ , 

{h.H] = ~{sgny)r-^{r ~ xf/^H = -^/2r-^/^ s\xv{4>/2)H . 
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8.2 Some other examples of integrability at fixed energy 

As discussed for example by Hietarinta [1], conformal transformations provide links between physically 
different systems which are integrable at some fixed energy. In particular if the original potential V is 
separable, V — Qi{X) + Q2{Y), then the transformed system has the potential 

V = \H\z)f[Qi{di{H{z)) + Q2mH{z)) - E] , (101) 

where z = x + iy and H{z) = X + iY. One of the results of the present work is that we have identified those 
conformal transformations for which the new potential in this situation is actually integrable at arbitrary 
energy (see table [IJ. Conversely, if the conformal transformation is not contained in table [T] then the 
resulting potential does not have a linear or quadratic invariant at arbitrary energies. Hietarinta considered 
conformal transformations of the forms H{z) — z™ (with m = —1, —2, ^,2), e^ and Inz. Note in particular 
that of these H{z) — z^" and H{z) — Inz produce systems which are integrable at arbitrary energy if the 
original potential is separable. In this subsection we give some further examples of simple systems which are 
integrable at a fixed energy. 

Consider first the polynomial function S{z) = iz^ . This is the simplest polynonfial which gives a potential 
which is not automatically integrable at arbitrary energy. The corresponding conformal transformation is 
given by w = ii{z) = 2~^/^(l — i) Inz or in terms of the real variables 

X = ^(0 + lnr), 

V (102) 

r=^(0-lnr). 

From the relation ((34| it then follows that the potential given by 

G = r-^[A{re'^) + Birer")] , (103) 

is integrable at zero energy for arbitrary functions A and B. For functions containing linear and quadratic 
terms the potential takes the form 

G = r-i (flie^ + aze-") + age^" + Uie'^'^ , (104) 

where the Oi (i = 1, . . . , 4) are arbitrary constants. 

As another example we take a function of the form S{z) = z^ where fc ^^ 0, 1, 2 is a real constant. This 
leads to w = H{z) = m^^z^"' where m = — fc/2 + 1 ^ 0, i, 1. The conformal transformation can then be 
written 

X — m^^r"^ cos(m0) , 

(105) 

— I 777, • / /i\ V ' 

= 171 r sin[mU) . 

The corresponding potential is 

G = r-''[A{X) + B{Y)] . (106) 

Choosing for example A{X) = aim^X'^ and B{Y) = a2m''F^ we have 

G = r-2+™(''+2) [ai cos^imO) + 02 sin" (to6I)] . (107) 

Specializing to the case to = 2 while keeping s arbitrary and using cos(20) = (a;^ — j/^)?-^^, sin(26') = 2xyr^^ 
yields finally the potential 

G = r^[ai{x^ - yy + a22'x'y'] , (108) 

which is therefore integrable at zero energy. 
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9 Concluding remarks 

It was shown in section |3] that confornial transformations given by analytic functions H{z) for which the 
condition 3{S"'(z)} = with S{z) — [iJ'(z)]^^ is satisfied give rise to the classical potentials which admit 
quadratic (or linear) second invariants at arbitrary energies. Conformal transformations which do not satisfy 
that condition give potentials which admit quadratic second invariants only at a fixed energy. Our approach 
unifies the description of quadratic invariants at arbitrary and fixed energies. In particular the integrability 
condition ([27)1 is valid for both types of invariants. It reduces to Darboux's classical condition for arbitrary 
energy invariants when Q'{5"'(z)} = 0. Whether a unification can also be achieved for third degree invariants 
or higher remains an open problem. 

An intriguing aspect of the integrability condition ()27p is the possibility that for a given potential function 
G = E — V there could exist a family of solutions S{z, E) with a continuous dependence on the energy. This 
would lead to new families of integrable potentials with energy dependent quadratic invariants. At this point 
we cannot exclude the existence of such solutions of the integrability condition. A related result was given 
by Hietarinta [3] who showed that the potential x/y is in fact integrable by energy dependent invariants 
which are certain transcendental functions of the momenta. 

For a given potential V the integrability condition ([?7|) with G — E ~ V can be used to determine 
energy values for which there exists an invariant of at most second degree. This involves solving a linear 
differential equation of the second order. For linear invariants it is sufficient to solve the linear equation 
([56)1 . We consider this possibility to test for linear and quadratic integrability at fixed energy to be an 
important application of the geometric approach to Hamiltonian dynamics. Another approach is to look 
for conditions involving curvature invariants such as the technique used in section 18.11 to find cases with 
additional invariants. In fact, conditions involving curvature invariants for (l+l)-dimensional models have 
recently been found (using two different approaches) for potentials which do not require the manifest linear 
invariant present in the Kepler problem [191 HQ] ■ There are indications that at least the approach of [19] 
can be generalized to incorporate quadratic invariants as well. It is of considerable interest to develop these 
techniques and use them to look for fixed energy invariants of physically interesting models such as the 
Henon-Heiles potential and others. 
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